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Abstract. We establish a new type of local asymptotic formula 
for the Green's function Qt{x, y) of a uniformly parabolic linear op- 
erator dt~L with non-constant coefficients using dilations and Tay- 
lor expansions at a point z = z(x,y), for a function z with bounded 
derivatives such that z(x, x) = x £ M. N . For z(x, y) — x, we recover 
the known, classical expansion obtained via pseudo-differential cal- 
culus. Our method is based on dilation at z, Dyson and Taylor 
series expansions, and the Baker-Campbell-Hausdorff commutator 
formula. Our procedure leads to an elementary, algorithmic con- 
struction of approximate solutions to parabolic equations which 
are accurate to arbitrary prescribed order in the short-time limit. 
We establish mapping properties and precise error estimates in the 
exponentially weighted, L p -typc Sobolev spaces W^' P (M. N ) that ap- 
pear in practice. 
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1. Introduction 

We establish a new type of local estimate for the Green's function 
of a uniformly parabolic linear operator with non-constant coefficients 
that do not depend on time. More precisely, we consider second-order 
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differential operators L of the form 

N N 

(1.1) Lu(x) := aij(x)did,ju(x) + bk(x)dku(x) + c(x)u(x), 
t,j=i fe=i 

where x = (xi, ...,Xjv) £ R"^, <9fc : = gf-> and the coefficients a^, 6j, and 
c and all their derivatives are assumed to be smooth and uniformly 
bounded. (We then write , bj,c e \M. N } and we denote the class 
of these operators by L.) We also assume that L is uniformly strongly 
elliptic, namely that there exists a constant 7 > such that 

N 

(1-2) E a ^)^>7ll£l| 2 , llef:=E^ 

ij i=l 

for all (£, x) € R^ x R-^. We define the matrix A(x) := [ay(x)], which, 
without loss of generality, we can assume to be symmetric. In view of 
the applications we are interested in, we take the coefficients of L to be 
real- valued. The set of operators L eh satisfying ( 11.21) will be denoted 
by L 7 . 

We study the short time asymptotic of the initial value problem (IVP) 
for the parabolic operator d t — L: 

d t u(t, x) — Lu(t, x) = g(t, x) in (0, 00) x 

u(0,z) = /(z), on {0}xM N , 



1.3) 



for u, f, and g in suitable function spaces. In view of Duhamel's 
principle, we may assume g = 0. 

When b(x),c(x) 7^ in (11.11) . the corresponding parabolic equation 
dtu — Lu = is collectively referred to as a Fokker-Planck equation. 
Fokker-Planck equations arise in many applications, for example in 
statistical mechanics [T5l [26] . and more generally in probability. 

We can also replace ~R N with a manifold of bounded geometry [181 
|4"T] . which thus allows us to treat also some degenerate elliptic operators 
L. In particular, the approach in this paper can be extended to the 
case of operators of the form d t — {ax 2 d^ + bxd x + c) acting on R f x R x+ 
and to other operators that appear in practice. This extension is work 
in progress [18]. See also below for a more detailed discussion of this 
point. 

It is known that there exists Q L 6 C°°((0, 00) x R N x R N ) such that 

(1.4) u(t,x)= [ g L (t,x,y)f(y)dy, t > 0, 

is a solution of the above equation, and it is unique if / and u satisfy 
certain growth conditions, specified later (see for instance [22], page 
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237). We will often write Q L (t,x,y) = Q^(x,y). In case we have 
uniqueness, we shall also use the notation u(t) = e tL f. The operator e tL 
is then called the solution operator of the problem (11.31) . and its kernel 
the Green's junction, or fundamental solution of L, or conditional 
probability density in applications to probability. 

For L with constant coefficients and for a few other cases, one can 
explicitly compute the kernel Q L . In general however, it is not known 
how to provide explicit formulas for Q L , though there is a large litera- 
ture on developing methods to obtaining good asymptotic formulas for 
the Green's function for t small and x close to y. For example, inter- 
preting the operator I as a Laplace-Beltrami operator on a manifold 
plus lower order terms, lead to formal asymptotic expansions of the 
form 

Gt (*, y) = (G (0) (x, y) + (*, y)t + ( X , y )r +...), 

as t — > + , where d(x, y) is the geodesic distance between x and y and 
Qv\x,y) are smooth functions in x and y. Among the vast literature 
we refer to [H SHI E31 S21 SHI [S3 [S3 EZ] , (see also [22 [23 SSI ED] for 
a pseudo-differential operator perspective). However, one difficulty in 
the practical implementation of this geometric approach is that, except 
again in special cases, there is no closed form solution to the geodesic 
equations used in defining d(x,y), which thus needs to be accurately 
approximated or computed numerically. 

A related short-time asymptotic approach uses oscillatory type inte- 
grals, which gives: 

(1.5) G L (t,x,y) ~ 52tW Pi (x,t-W{x-y)) e - ^§^ , 
i>o 

as t — > 0+, where Pj(x, w) is a polynomial of degree j in w, and A(x) : = 
[aij(x)]. (We follow here Taylor [50l Chapter 7, Section 13], where an 
asymptotic parametrix for the heat equation on compact manifolds 
was constructed.) Finally, we mention the recent approaches in [1] 
using multivariate Hermite expansions, and in [TD] using an alternate 
construction of a parametrix approximation. 

In our paper, we devise a new, elementary method to obtain asymp- 
totic expansions similar and even more general than (11.51) . Our method 
is based on dilating the coefficients of L around a point z with ratio 
t 1//2 , then expanding in a Taylor series in t. We regard this expansion 
as a perturbation of the operator Lq obtained from L by freezing co- 
efficients at z. We then use a Dyson-series perturbative expansion to 
approximate the heat-kernel of L. The Dyson-series expansion turns 
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out to be explicitly computable using the Baker- Campbell-Hausdorff 
commutator formula. We call the resulting method the Dyson-Taylor 
commutator method. We think that our method is more accurate and 
more stable in practical implementations pJH [18] . 

The main goal is to provide an explicit, algorithmic method to com- 
pute each term in the expansion, while at the same time obtain sharp 
error bounds in both weighted and unweighted Sobolev spaces. We do 
not work on compact manifolds, rather in R N , so that the error needs 
to be globally controlled (see below for a connection with operators 
on non-compact manifolds of bounded geometry). In particular, our 
approximation is valid uniformly in x and y, provided t is small enough. 

We think that our method, the Dyson- Taylor commutator method, 
is more elementary than the ones found in the literature and since it 
relies on an iterative time-ordered perturbative formula for the solution 
operator e tL , Equation (12.151) . a parabolic rescaling argument, and a 
suitable Taylor's expansion of the coefficients of L (equation l3.8l) . Since 
the iterative formula is obtained via repeated applications of Duhamel's 
principle, we could also treat certain classes of semilinear equations 
following Kato's method, which allows to take rougher data as well 
(see [3S] in the context of the Navier-Stokes equations). We remark 
here that a similar parabolic scaling combined with Taylor expansions 
has been used in obtaining a short-time expansion for stochastic flows 

(see nnna). 

Our main result is the following theorem. We introduce the weight 
(x) = (1 + \x\ 2 ) 1 / 2 . Below, W™' p := W™' P (R N ) is the exponentially 
weighted Sobolev space defined by 

W™' P (R N ) := {u : R N — > C, d^(e a{x) u(-)) G L P (R N ), \a\ < m}, 

for 1 < p < oo, m G Z + , and a G R. (See also Equation (12.41) ). When 
a = 0, we recover the usual Sobolev spaces. The need to consider 
exponentially weighted spaces arises in applications to probability, in 
particular in stochastic volatility models. For instance, after making 
the substitution x = e y , the payoff usually associated with the Black- 
Scholes equation (equation (11.81) below) belongs to W™' p with m = 1, 
a < — 1, and p large. We also denote 

(1.6) G(z; x) = (4ir)-W det(A(z))- 1/2 e- xTAiz) ~ lx/ \ 

where z is a given point in R . It is interesting to mention that the 
Black-Scholes equation fits into the framework of manifolds with cylin- 
drical ends, to which the results of Krainer [37] apply. Manifolds with 
cylindrical ends are the simplest examples of manifolds with bounded 
geometry. 

4 



To approximate the value of the Green function Q t (x,y) at some 
point (x,y), we will use a Taylor-type expansion at the point z of a 
suitable parabolic rescaling of the coefficients of L, which, however will 
be chosen depending of x and y, z = z(x,y). Typically z(x,y) = 
Xx + (1 — X)y, for some fixed A, but we can allow more general choices. 
Namely, we shall say that z(x,y) is admissible if z(x, x) = x and all 
derivatives d a z are bounded for a^O. 

Theorem 1.1. Let fi £ Z+, L £ L 7 , z = z(x,y) be an admissible 
function. Let ty e (z,x,y) = a a ,p(z)(x — z) a (x — y) 13 , \a\ < I, (5 < 
3^, a a ,i3 £ C^°(1R 7V ) ; be the functions provided by the Dyson-Taylor 
commutator method explained in the second half of this Introduction. 
Define for each integer < £ < fi, 

g^\x,y) :=t-^jrt^^Uz+ 

where z = z(x,y). Define the error term S^'^ in the approximation of 
the Green's function by: 

e tL f{x) = [ g [ r\x,y)f(y)dy + & +1 ^S [ r ] f(x). 
Jr n 

Then, for any f £ W™' P (R N ), a £ R, m > 0, 1 < p < oo, we have 

(i-7) \\£^ z] f\\ wr+k , P < ct-wwfWw?*, 

for any t £ [0,T] ; < T < oo, k £ Z +; m£/i C independent of 

te[o,T]. 

The function Q^ J, ' z \x,y) will be called the fith-order approximation 
kernel for the solution operator e tL . 

See Subsection [TTTJ at the end of this introduction for a more detailed 
description of how the approximation kernel q\^' z \x, y) is obtained. In 
Section H] we give an explicit iterative construction of the functions < $ e . 
The main interest is, of course, in the derivation of the approximation 
kernel. However, without good error estimates, this kernel will not be 
of great use in practice. 

For the particular choice z(x, y) = x, we have checked that the first 
few polynomials pj(x,x — y) := ty 3 ! (x, x, y) coincide with the ones in 
the expansion (jl.5p above (see [501 Chapter 7, Section 13]). Our result 
is more general, however. We discuss in [18] different choices of the 
additional function z = z(x,y) in the framework of the usual Black- 
Scholes equation. It turns out that the choice z = x is not always 
the most appropriate. In fact, for the Black-Scholes equation and fi = 
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0, the choice z = Jxy can lead to a better approximation, whereas, 
surprisingly, the choice z = (x+y)/2 yields worse numerical results than 
simply choosing z — x. This addional accuracy obtained for a suitable 
choice of z(x,y) ^ x and in view also of the mid-point quadrature 
rule (which leads to a higher order of convergence) justifies the extra 
generality of including arbitrary admissible z in our method. 

A localization procedure as in [?T] will allow us to pass from opera- 
tors on Mr to operators on manifolds M of bounded geometry (again 
following [T7] and [2]). More precisely, our results will extend to op- 
erators of the form L = a ijdidj + bidi + c, defined on a subset 
Q of M. N such that the coefficients are bounded in normal coordinates 
with respect to the metric g = cf-'dxidxj, assumed to be com- 
plete of bounded geometry on Q. Here the matrix [a y ] is the inverse 
of the matrix A, following the usual convention. Such metrics arise 
naturally when resolving boundary singularities, (see [13] for a sys- 
tematic treatment of heat calculus on manifolds with boundary). We 
refer to [2j [T7J [351 E3] for recent papers dealing with partial differen- 
tial equations on manifolds with metrics of this form. In particular, 
we can deal with certain operators having polynomial coefficients such 
as those arising in probability and its applications, for example in the 
Black-Scholes option pricing equation [9] 

(1.8) Lu(x) = ax 2 d 2 ,u(x)/2 + r(xd x u(x) — u(x)), 

where in this context t is the time to option expiry. Our results also 
apply to differential operators arising in stochastic volatility models 
(c.f. [3J El ESI [27] ED ESI SO]). On the other hand, our results apply to 
operators of the form x 2f3 d 2 , < f3 < 1 only locally. A good framework 
for obtaining differential operators with unbounded coefficients that 
satisfy our assumptions is that of Lie manifolds [2J . This point will be 
discussed in detail in [18] . Explicit calculations and concrete, practical 
applications of our method will be given in [191 [20] • 

In addition, our methods generalize to operators with time-depdendent 
coefficients, satisfying certain conditions. This extension is addressed 
in a forthcoming paper |21j . 

We conclude this first part of the introduction with an outline of 
the paper. In Section [2J we define the weighted and regular Sobolev 
spaces of initial data for the parabolic equation and introduce the class 
of operators L under study. We also briefly discuss mapping properties 
of the semigroup generated by L and use them to justify the Dyson (or 
time-ordered) perturbation expansion of e tL . In Section [3j we exploit 
local in space and time dilations of the Green's function together with a 
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certain Taylor expansion of the operator L to rewrite the perturbation 
expansion as a formal power series in s = \Jt. In Section HI we employ 
commutator estimates to derive computable formulas for each term in 
the expansion. This leads to the Dyson- Taylor commutator method to 
determine the functions ^ used in Theorem 11.11 

Finally, in Section [5J we rigorously justify our expansion and derive 
error bounds in time by means of pseudodifferential calculus. 

Acknowledgments We thank Andrew Lesniewski and Michael Taylor 
for sending us their papers and for useful discussions. We also thank 
Richard Melrose for carefully reading our paper and for pointing out a 
possible missunderstanding in an earlier version of this paper. 

1.1. The approximate Green function. We close this Introduction 
by describing in more detail the Dyson- Taylor commutator methodused 
to define approximation kernel Q\ n ' z \ Given an operator T with smooth 
kernel, we denote its kernel by T(x,y), as customary. 

We consider a uniformly strongly elliptic differential operator L of the 
form L := £)ij=i a^didj + Y^k=i b kdk + c, where a ih b k ,c G C£°(R N ) 
are real valued. Given a fixed point z G WL N and s > 0, we define 
L s,z ■= Y^i,j=i a ij Z ( x )didj + s Y^i=iK' Z ( x )^i + s 2 c s > z (x), where for a 
generic function / we set f s ' z (x) = f(z + s(x — z)). Hence, z acts 
as a fixed dilation center. We then Taylor expand this operator in s 
around to order n: 

n n+1 

(1-9) L s > z = £ s m L z m + V£ x = * m L Z m, 

m=0 m=0 

where Lj, < j < n, are differential operators with polynomial coeffi- 
cients that do not depend on s, whereas L z n+l has smooth coefficients 
that, however, do depend on s (although this dependence is not shown 
in the notation). Hence, V^ii = s n+l L z n+1 is the remainder of the 
Taylor expansion. The order n will be chosen later. In particular, we 
observe that 

(1.10) Ll = ^a ij (z)d i d j . 

hi 

For any fixed, positive integers k < n + 1 and £, we shall denote 
by 2lfc/ the set of multi-indexes a = (cki, a.i, ■ ■ ■ , &k) £ such that 
\a\ := Yl a j = (Hence, k < £.) Then, for each multi-index a = 
(ai, . . . , ttfc) £ 2lfc/, k < n, we introduce 

(1.11) A a , z := [ e T ^L z ai e T ^L z a2 ---L z ak e T ^dr. 
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Observe that k, £ are unique given a. 

The main point is that the operators A a can be computed explicitly 
as follows. Let us denote by T>(a, b) the vector space of all differen- 
tiations of polynomial degree at most a and order at most b. (By 
polynomial degree of a differentiation we mean the highest power of the 
polynomials appearing as coefficients.) Then, for any Lq G X>(0, 2) 
that is uniformly strongly elliptic and for any L m G D(m, 2), we 
have a differential operator P m (L , L m ;6,x,d) given by the formula 
e eL "L m = P m (L ,L m ;6,x,d)e eL ° 1 where 9 > (see Lemma 03]). Let 
Efc be the unit k- dimensional simplex. Next, for any given mult i- index 
a G 2l M with k < n, we define V a (x, z, d) := [jLi p ^( l o> 1 - 
(Ji,x,d)da. Then 

(1.12) A QjZ = P Q (x, 2; ,^)e i o, 

where the product is the composition of operators and V a is a differen- 
tial operator of order 2k + £ in x and polynomial degree < i in (x — z) 
(see Lemma [4. 6 1) . 

Since z is arbitrary, but fixed at this stage, if Lq is the operator in 
f fLTOj) then e L o(x,y) can be explicitly calculated and it agrees with the 
function G(z, x — y) introduced in equation (jl.6p . Therefore, it can be 
easily seen from (11.121) . that 

A Q , z (x, y) = ^(z, x, y)G(z; x-y), 

for some 9($ (z, x, y) = J2 a u,p(. z )(. x ~ z) a (x — y) 13 , \a\ < £, [3 < 3£, 
a Q ,/3 £ C%°(M. N ). In particular, all A Q ^ are operators with smooth 
kernels, thus denoted A a ^ z (x, y). 

We will show that e tL as well is an operator with smooth kernels, 
henceforth denoted Q^(x,y). 

Let us fix for the time being a smooth function z : M. N x M. N —>■ ~R N 
whose properties will be made precise below. (Two typical examples are 
z(x, y) = (x+y)/2 and z(x, y) = x, which suffice in many applications.) 

Our approximation will be obtained by combining Lemma [3.21 with 
the perturbative estimate of Equation (13.181) at some point z = z(x, y) 
using the dilation with center z(x,y) and denoting s 2 = t. Then, for 
any /i > n, we define 

(1.13) Q [ r\x,y) := s- N e L *(z + s- 1 (x-z),z + s-\y-z)) 

+ Yl s e A a , z (z + s- 1 (x-z),z + s- 1 (y-z)). 

£=1 k=l ae2t fc , £ . 
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The operator L is not shown explicitly in the notation Q t (x, y), al- 
though Qf' ,z \x,y) does depend on L. This is not likely to cause any 
confusion, since L is usually fixed in our discussions. 

The justification of the above definition for the approximation is that 
a Dyson series expansion of order n + 1 gives us 

(1.14) g t (x,y) :=s- N e L »(z + s- 1 (x-z),z + s- 1 (y-z)) 

(n.+l) 2 max{£,n+l} 

+ Y Y Y s l K^ z {z + s~\x - z),z + s^ l {y - z)), 

1=1 k=l ae%u,i 

where for k — n + 1 and a = (a±, . . . , a^) G 21^, we introduce 

(1.15) A a , z := [ e^L z ai e T ^L z a2 ---L z ak e^ z dr. 

The difference between equations (11.131) and (11.141) is that the sum in 
the first equation contains exactly the terms with s , £ < n, from the 
second equation. The difference between equations (11.111) and (11.151) 
is in the last exponential. Note that A Q)Z does not depend on s if i = 
\a\ < n, but it may depend on s otherwise. In any case, all the terms 
Aq,^ that depend on s will be included in the error term. All terms A Qi2 
with fi < \a\ < n, which do not depend on s, will also be included in 
the error. We remark that the error term is never computed explicitly, 
as only the /zth order approximation kernel is needed. Therefore, while 
H will usually be small in applications, we can take n as large as needed 
to justify the error bounds of Theorem ll.il In Section [51 we will show 
that n > fi + N — 1 suffices. 

2. Preliminaries 

We begin by discussing in more details the class of second-order 
operators L of the form (11.11) that are the focus of our work. Below we 
set 

(2.1) C°°(R N ) := {/ : R N — > C, d a f bounded for all a }. 

Definition 2.1. We shall denote by L the set of differential operators 
L of the form 

N N 

(2.2) L: =Y ai & d i + Y bkdk + c ' 

i,j=l k=l 

where aij,bk,c G C£°(M. N ) are real valued. We shall denote by L 7 
the subset of operators L G L satisfying the uniform strong ellipticity 
estimate (11.21) with the ellipticity constant 7. We let A = [a^] and 
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assume additionally that A is symmetric, which can be achieved simply 
by replacing A with its symmetric part, since this does not change our 
differential operator. 

The above definition can be extended to operators on manifolds of 
bounded geometry M (see [T71 |4"TI H7j). For example, when M = R N 
with the Euclidean metric, the class L considered in [JT] coincides with 
the class L considered in this paper. 

In what follows, we denote the inner product on L 2 (R N ) by [u, v ) = 
J RN u(x)v (x)dx. Let us denote (£) := (1 + l^] 2 ) 1 / 2 and let u be the 
Fourier Transform of u. We also recall the definition of and some basic 
facts about L p -based Sobolev spaces W r,p (R N ) . For 1 < p < oo, r G R: 

(2.3) W r ' p (R N ) := {u : R N C , (£) r u G L P (R N )} 

= W r,p (R N ) := {u : R N -»• C , (1 - A) r/2 u G L P (R N )}, 

If r G Z+, 

W r ' p {R N ) = {u : R N -> C, <9 a u G L P (R N ), \a\ < r}. 

Since the dimension N is fixed throughout the paper, we will usually 
write W r > p for W r ' p {R N ). When 1 < p < oo, the dual of jy r ' p is the 
Sobolev space W~ r ' p with l/p+l/p' = 1. 

We are interested in considering the initial value problem (II. 3p in the 
largest-possible space of initial data / that includes the typical initial 
conditions that arise in applications and where uniqueness holds. We 
therefore introduce exponentially weighted Sobolev spaces. Given a 
fixed point y G R N , we set (x) w := (x — w) = (1 + \x — w] 2 ) 1 ^ 2 and 
define W™*{R N ) for m G Z+, a G R, 1 < p < oo, by 

(2.4) W^K") := e" a<a;>1 "W' p (R JV ) 

= {u : R N -> C, ^(e a ^-w(-)), G //(R^), |a| < r}, if r G Z+, 
with norm 

a|<m 

When it is clear from the context, we may drop the subscript w from the 
above notation. We observe that W™' p = W™£ = W m,p . The spaces 

W£' p and Wla'w are naturally duals to each other if l/p+l/p' = 1. 

A crucial observation is that, for any L G L 7 and any a G R, the 
operators L a := e al ^ w Le~ a ^ w are also in L 7 . They moreover define 
a bounded family in L 7 if a is in a bounded set, while w is arbitrary. 
Since proving a result for L acting between weighted Sobolev spaces 
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* s the same as proving the corresponding result for L x acting 
between the Sobolev spaces W s,p = W^, we may assume that a = 
and w is arbitrary. In particular, L : W^ 2 ' p — > is well defined 

and continuous for any a and w, since this is true for a = 0. 

In fact, it will be crucial for us to establish mapping properties that 
are independent of w . This will be the case in all estimates below, unless 
stated otherwise. One of the most important example is provided by 
Corollary 13.91 Moreover, the spaces W™* do not depend on the choice 
of the point y (although their norm obviously does). Because of this 
observation, we shall often omit the point w from the notation, when 
this does not affect the clarity of the presentation. 

We begin by recalling some properties of L and the associated solu- 
tion operator e tL to the initial value problem (11 .3j) . 

2.1. Mapping properties. Given a Banach X and an interval I of the 
real line, we shall denote by C(I,X) the space of continuous functions 
u : I — > X. By C k (I,X) we shall denote the space of functions u G 
C(I,X) such that u {j) G C(I,X) for all < j < k. We assume that 
X C Ll oc (M. N ), and that L is a closed unbounded operator on X with 
domain V(L) C X. 

Let g G C([0, oo),X). By a classical solution in X of (11.31) we mean 
a function 



(2.5) u G C([0,oo),X) nC 1 ((0,oo),X)nC((0,oo),P(L)), 



such that d t u(t) = Lu(t) + g(t) in X for all t > and it(0) = / in X. 
(The domain of L is given the graph norm |||it||| := \\u\\ + which 
makes T>(L) a complete normed space, since we have assumed that L 
is closed and X is complete.) In particular, tt(0) = / must belong to 
the closure of T>(L) in X. In the case of interest here, if X = W^ ,p , 
then V(L) = W^ +2,p , which is dense in X. 

In view of Duhamel's formula (which will be justified below), we 
can assume g = in Equation (11.31) . We shall take our Banach space 
where the solution is defined to be X = L p for some arbitrary, but 
fixed, p G (1, oo) and a > 0. Then Equation (II. 3p becomes 



Let us notice that if / G C°°(R Ar ) also, then we recover Equation (II. 3p . 
The growth condition u(t) G L p is needed, however, in order to insure 
uniqueness. 

A family of (bounded) linear operators U(t) on X, t > 0, will be 
called a C° or strongly continuous semigroups of operators if U(t) forms 



(2.6) 




in L p (R N ), 

f e L p (R N ). 
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a semigroup in t and U(t)u — > u in X as t — > 0+. This last property 
shows that the function [0, oo) 9 t — > U(t)f G X is continuous for any 

/ex. 

We shall need the following standard result. Recall the subset L 7 C 
L introduced in Definition 12.11 

Lemma 2.2. (i) Let L G L 7 , i/ien t/iere exists a constant C > suc/i 

7(Vn, Vm) - C(u, u) < -{Lu, u) < C(Vu, Vu) + C(u, w). 

(^'^ T7ie norm |||f|||2m := IMIlp + ll-^ m ^IUp i> s equivalent to the norm 
| • ||jy2m,p on W 2m ' p (R N ) , for any m G Z + and 1 < p < oo. 

Proof. (Sketch.) (i) follows from a direct calculation. See [17], [5Tjor 
gl] for (ii). □ 

It follows from this lemma that L : iy 2,p — > LP is a closed, densely 
defined unbounded operator on LP . This technical fact is important 
because it is often needed for the general results that we will use below. 

For the sake of clarity and completeness, we include here a quick re- 
view and some proofs of the main properties of the semigroup generated 
by L. Our proofs also serve the purpose of justifying the perturbative 
expansion described in Section 12.21 which is discussed extensively in 
the literature, but usually not in the setting that we need. Further 
details can be found in [3H [3U 05]. Below p G (1, oo) and 7 > will 
be arbitrary but fixed, and the constants appearing in the estimates 
depend on p and 7, but not on L G L 7 . 

Proposition 2.3. Let 06K, 1 < p < 00, and L G L 7 . 

(i) For each f G W% ,p , the problem (11.31) has a unique classical solu- 
tion 

u g c([o, 00), Li) n ^((0, 00), Li) n c((o, oo), w^ p ). 

(ii) Let e tL f := u(t), then we have e tL W^ ,p C W^ ,p and, moreover, 
\\ etL f\\wa' p — Ce ut \\f\\wa ,p > f or a constant C independent of r, a, 
and L G L 7 in bounded sets.. 

Proof. We can assume a = 0, as explained above. Lemma I2T21 (i) gives 
that L satisfies the assumptions of the Hille-Yosida theorem [23l 1381 EE5] . 
and hence e tL is defined, is a C° semigroup, and u(t) := e tL f is indeed 
a classical solution. This proves (i). 

It also follows from standard properties of C°-semigroups in Banach 
spaces that ||e tL /||i P < Ce^ for some constants C > andcj G K 
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independent of L G L 7 . To prove (ii), we then notice that 



(2.7) \\e tL f\\ w ^ P < C\\\e tL f\\\ 2m , p = C(\\e tL f\\ LP + \\L m e tL f\\ LP ) 
= C{\\e tL f\\ LP + \\e tL L m f\\ LP ) < Ce^{\\f\\ LP + \\L m f\\ LP ) 

= Ce^\\\f\\\ 2m <Ce^\\f\\ w ^ P , 

with constants depending on m, p, and L, but not on t. Though L 
may not be self-adjoint, the adjoint L* is an operator of the same type, 
in the sense that L* G L 7 . Hence the estimate above holds for L*, 
with possibly different constants. We can then extend Equation ( 12.71) 
to W~ 2m ' p , m G Z+, by duality and to any W r,p by interpolation (see 
for example [8j [52] for results on interpolation). This completes the 
proof. □ 

From now on we shall denote by e tL the C°-semigroup generated by L 
on L p = e "^ z L p (R N ) , with p and a determined by the context (usually 
arbitrary, but fixed). 

We recall that for / G V(L), the map t -> e tL f is in ^([0,00), X) 
and dte tL f = e tL Lf = Le tL f. For any two normed spaces X and Y, we 
denote by B(X, Y) the normed space of continuous, linear operators 
T : X — > y with norm ||r||x-»y- When X = Y", we shall also write 
||T||x := ||T||x^x and B(X) := <B(X, X). The identity operator of any 
space will be denoted by 1. 

Lemma 2.4. Let L G L 7 . M^e /icrae ||e tL — 1|| ^+2,5^^,, < Ct, for 
any t G (0,1]. In particular, [0, 00) 3 t -> e* L G° S(w/ +2 >p W"*'*') zs 
continuous. 

Proof. We have e* L / — f = f^e sL Lfds for any / G W 7 ^, by standard 
properties of C°-semigroups. Lemma [2.31 ((n]) then gives 

\\e tL f - f\\ w s, < I \\e sL \\ w , P \\Lf\\ w , P ds < Cte^\\f\\ w s +2 , P , 
Jo 

which proves the first part of the result. 
Let now t\ >t 2 . Then 

\\„tiL „t 2 L\\ ^ \\Jti-t 2 )L ill lUfeill 

||e -e \\ w s+2, P ^ w s, P < ||e v - lH^a+^^.p ||e || w .,p. 

This completes the second part of the proof. □ 

Remark 2.5. Let 5 G (0,2]. Then an interpolation argument gives 
\\e tL — l\\ wS +s,p_^ w s,p < Ct s / 2 , for any t G (0,1]. Hence the function 
[0, 00) 3 t A e tL eB(W° +s ' p , W^ p ) is also continuous. 
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We discuss smoothing properties of e tL , it is convenient to first as- 
sume L* = L, that is that L is self-adjoint. This will require us to set 
a = in our weighted Sobolev spaces W* ,p = W^ ,P (M N ). This assump- 
tion will be removed later on. The following result is known, we sketch 
a proof for completeness. (See for example [15] and [H] in the more 
general case of manifolds with bounded geometry.) 

Corollary 2.6. Let t > 0. There exist constants C ryS > such that, 
for any LeL 7 with L = L* : 

(i) \\e tL f\\w-p(m) < C r ^ s - r)/2 \\f\\ W s, P{R N), r >s real. 

(ii) There exists GH x iV) e C°°((0,oo) x R N x R N ) such that 

(2.8) e tL f(x)= [ G}(x,y)f(y)dy. 

Jr n 



Proof. The part can be proved using resolvent estimates and a 
scaling-in-time argument. Part (JD} folllows from the Schwartz ker- 
nel theorem (see for example [501 Chapter 7]), since from e tL maps 
compactly supported distributions in £'(W N ) to smooth functions in 
C°°(M. N ). In fact, if we denote by <, > the duality pairing between C°° 
and £', we explicitly have: 

(2.9) Gt , (?,v)=<6 x ,e tL 6 v >, 

where 8 Z , z e Mr, represents the Dirac delta distributions supported 
at *(i.e., &(/) = /(*)). □ 

We now proceed to eliminate the assumption that L* = L in the 
above result. First, let us notice that if L, Lq e L 7 , and if we denote 
V = L — L and g(t, x) = Vu(x, t), then (11.31) becomes 



(2.10) 



d t u — Lqu = g in (0, oo) x M. N , 

u(0,x) = f(x) on {0} x 1^. 



It is well-know that applying Duhamel's formula, gives a Volterra in- 
tegral equation of the first kind for u. If L = Lq, the solution of the 
integral equation is a classical solution of (12.101) . in fact, it is enough 
that e tL ° generates an analytic semigroup (see [151 Theorem 2.4, page 
107]). For simplicity, we want to avoid using the theory of analytic 
semigroups, and rather use instead Corrolary 12. 61 

Lemma 2.7. Let us assume that g G C([0, oo), L p ). Then the classical 
LP -solution of the problem (I2.10p is given by 

(2.11) u (t) = e tLo f+[ e {t ' T)La g(T)dT. 

Jo 
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Assume that L £ L 7 , and let L = (L* + L)/2. Then, the classical 
LP -solution u(t) =: e tL f to the problem (12.61) is given by: 

(2.12) e tL f = e tL °f+ [ e {t ~ T)Lo (L - L )e TL fdr, 

Jo 

for any initial data f £ L p , 1 < p < oo. 

Proof. Let us notice that it is denned since e^~ T ^ L °g(T) is continuous 
in r. For / £ W 2,p and g £ C([0,oo),W 2,p ), the function it is also 
differentiable and 

u'(t) = L Q e tL °f+ [ L e^ L "g(r)dr + g(t) = L Q u(t) + g(t). 
Jo 

Since C([0, oo), W 2 ' p ) is dense in C([0, oo), L p ) (for the topology of uni- 
form convergence on compacta), this proves the first part. 

To prove the second part, let us chose / £ W 1,p , then the function 
[0, oo) 3 t — > (L — L )e rL f = g{r) £ LP is continuous (since L — L is 
in L and has order at most one), and hence we can apply the first part 
to obtain the formula (I2.12p . 

In general, using Proposition 12.31 part ([n]), and the fact that V £ L 
and has order at most one, we obtain by Corollary 12.61 P ar t © that 
|| e (i-r)Lo( L _ L )e TL f\\ LP < Ct- l l 2 \\f\\ LV , so that the integral on the 
right hand side of (12 . 1 2[) is defined and continuous in / £ L p . Since 
the left hand side of (12.121) is also continuous in / £ L p , the result then 
follows by continuity and by density of W 1,p in LP . □ 

Remark 2.8. The integral J * e^~ T ^ L °(L — L^)e rL fdr is defined either 
as a Bochner integral (for the definition of the Bochner integral see 
e.g. [44]) or as the limit lim^o e e^~ T ^ L ° (L — L )e TL fdr of Riemann 
integral for continuous functions. 

We now extend Corollary 12.61 to non self-adjoint operators L and to 
the exponentially weighted spaces W*' p . 

Proposition 2.9. Let L £ L 7 arbitrary. We have e tL W^' p C W^' p for 
all r, s,a £ M, 1 < p < oo, and t > 0. Let r > s, then 

\\ U (t)\\ W r, P <Ct^ 2 \\f\\ W , P , t G (0,1]. 

The constant C above is independent r, s, a, p, and L, as long as they 
belong to bounded sets. 

We recall that W^' p is independent of the choice of the point z (see 
Equation (12.41) ). The constant C in the above proposition is also inde- 
pendent of z £ M N since the family e a ^ z Le~ a ^ z is uniformly bounded 
in L 7 for z £ M N and a in a bounded set. For this reason, we shall 
sometimes drop the index z from the notation (x) z . 
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Proof. As discussed above, we may assume that a = 0. Also, note that 
we already know that e tL W^ p C W^ p for all r < s, so let us concentrate 
on the non-trivial case r > s. Let L := (L + L*) /2 and V = L — L . 
Then V G L is a differential operator of order at most one. By Lemma 

EH 



e tL = e tL + f e (t~r)L y e rL dT _ 

Jo 



Let us assume also that s < r < s + 1. Using also Proposition 12. 3[ 
part HU we obtain that the norm ||e tL ||iy s >p^w r >p of e tL as linear map 
W s ' p — > W r ' p can be bounded as 

|| C || \ya,p_ >\yr,p ^ || e ||^s,p_ >Wr,p 

+ /" ||e^* || yys-i,p^yyr, P || V || vi/ s 'P^iy s_1 'P II c T ^ 1 1| w s 'pd>T 
t 



-«)/2 



< C ^(*-r)/2 + y (t _ r) ( S -l-r)/2 dT ^ = ^-0/2(^1/2) < ^(r 

where in the last inequality we have used that < t < 1, and where 
C > is a generic constant, different at each appearance. The general 
case follows from this one as follows. Let 5 = (r — s)/m, for m > 
r — s. We first notice that ||e* i / m || H / s +o-i)<5,p^yi/s+j<5,p 
by the result that we have just proved, since 5 < 1. We then write 
we use the submultiplicative property of the norm 
to obtain He^il^^r* < C '{t / Wp^' = Ct^' 2 . □ 

In particular, Proposition 12.91 gives the existence of the Green's func- 
tion Q^(x,y) for any L G L 7 , defined again via formula (I2.9p . In par- 
ticular for t > 0, this kernel is a smooth function of x and y. We will 
also use the notation Q^(x, y) = e tL (x, y). The following corollary is a 
consequence of Proposition 12.91 

Corollary 2.10. Let L G L 7 and s,r G R be arbitrary. We then have 
that the map 

(0, oo) 3 t - e tL G W™) 
is infinitely many times differentiable. 

Proof. We have d k e tL = e tL L k , so it is enough to show that the map 
(0, oo) 3 t — > e tL G B(W*~ k,p , Wl' p ) is continuous. Now, for each 
5 > 0, let t > 5 > 0. Then e 5L maps W°~ k ' p to W^ +2 ' p continuously, by 
Proposition 12.91 Writing e tL = e^~ 5 ^ L e 5L and using the continuity of 
[5,oo) 3 t -> e^ 5 ^ G i3(iy a r+2 ' ?, ,iy a r ' ? '), by Lemma E31 we obtain the 
result. □ 
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See [131 [m HH] for more continuity properties of the semigroups gen- 
erated by second order differential operators (Schrodinger semigroups). 

Let us notice for further reference that for constant coefficient oper- 
ators, the Green's function can be determined explicitly. 

Remark 2.11. If L (11. ip is a constant coefficient operator 



(2.13) i°=E fl ^ + E^ 

i,j=l k=l 



c° 



and A := (a° ) is the matrix of highest order coefficients, assumed to 

satisfy a°- = a° i; we have the explicit formula 

(2.14) 

„ T 0, . tT , , e Ct ( 3! +6 t-i/) t (A )- 1 (x+b t-i/) 

* v /(47rt)«det(A°) 

2.2. Perturbative expansion. The purpose of this section is to ob- 
tain a time-ordered perturbative expansion of e tL , L e L 7 , in terms 
of e* L ° for a fixed element Lo G L 7 . Later, Lo will be obtained by 
freezing the highest-order coefficients of L at a given point 2 and drop- 
ping the lower-order terms. This expansion is the well-known Dyson 
series f2E \ I3"2"l I3"4] . Here, we concentrate on justifying this expansion in 
our setting and in obtaining global error estimates in weighted Sobolev 
spaces. 

For each k G Z+, we denote by 

S fc := {t = (to, 7i, ... , 7fc) € K fc+1 , 7} > 0, ^ Tj = 1} 

~ {a = (o"i, . . . , a k ) e R k , 1 > 01 > (7 2 > . . . (Tfe-i > <Tfc > 0} 

the standard unit simplex of dimension fc. The identification above 
is given by Uj = Tj + Tj+i + . . . + Tk- Using this bijection, for any 
operator- valued function / of ~R N we can write 

I f(r)dr= / / .../ f(l-ai,a 1 -a 2 , ■ ■ ■ ,a k -i-a k ,a k )dak ■ ■ .dai 

JE fe Jo JO Jo 

= / /(I — 01)01 — <72, • • • ,0fc_i — <7fc,CTfc)d£7 

Js fc 

We recall that, if g : [a, 6] — > X is a continuous function to a Banach 
space X, J g(t)dt is defined as a Riemann integral. We begin with a 
preliminary lemma. We further recall that B(X, Y) the Banach space 
of continuous, linear maps between two Banach spaces X and Y. 
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Lemma 2.12. Let Lj £ L 7 and let Vj be such that e b i^>Vj £ L, 
j = 1, . . . , k, for some b = (b\, . . . , bk) £ k £ Z + . Tnen 

$(r) = e Toi °yie riil . . . e T ^ L ^V k e T * L * , r £ S fc 

de/mes a continuous function $ : S fe -> /B(W*' P (R JV ), W^f.y^)) /or 
any a £ M and 1 < p < oo. 

Above we use the standard multi index notation |6| = X]j=i 

Proof. It is enough to prove that <3> is continuous on each of the sets 
V? := i T j > V(^ + 2)}, j = 0,...,k, since they cover S^. Let us 
assume that j = 0, for the simplicity of notation. 

By assumption and by Lemma 12.41 each of the functions 

[o, oo) 3 Tj _> Vj e^ £ B(w^, ii ;:; ';, ). 1 < j < *, 

is continuous. For a suitable choice of Cj and rj (more precisely, Cj = 
cj + i — bj + i, Ck = a, Tj = Tj + \ — 4, rfc = s), we obtain that the map 

[0,oo) fe 3 ( Tj ) =: r' - tf(r') := V^ lLl ...^ Lfe £ B(W a s ' p , W^tf*) 

is continuous. 

Corollary EIUl gives that the map r -> e r ° L ° £ i3(W a s " 4 6 f' p , W a r f |b| ) is 
continuous for t > l/(fc + 2). This proves the continuity of $ on Vo 
and completes the proof of the lemma. □ 



By iterating Duhamel's formula in Lemma \2.7\ we obtain a time- 
ordered expansion of e tL . 

Proposition 2.13. Let d £ Z+. Then, for each L, Lq £ L 7; 

e tL = e tL 0+t I e tT L y e tTiLo dr 



(2.15) 



+ t 2 / e^oioyginioygtraLo^ + . . . + 
+ t d / e tT Loy e tTiL _ _ _ g^-i^oygir d L orfr 

+ / e tToLo Ve tTlLo . . .e tTdL °Ve Td+lL dT, 

where V = L — Lq, and each integral is a well-defined Riemann integral 
of a Banach valued function. 

The positive integer d will be called the iteration level of the approx- 
imation. As d — * oo, formula (12. 15[) above gives rise to an asymptotic 
series (Dyson series, see [281 E2l [M] and the references therein). 
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Later in the paper, V will be replaced by a Taylor approximation of 
L, so that V will have polynomial coefficients in x, so we have included 
this case in the lemma above. 

Proof. Recall that Lemma [2.71 gives 



e tL _ e tL 



f e^OLoy e (L dc = r e t(l-T)L y e trL td ^ 

Jo Jo 



with the substitution £ = tr. This is in fact our result for k = 1. 

The result for any p then follows by induction using the above for- 
mula. 

Recall that on each simplex E p , we denoted <Tfc = + rfc+i + . . . + r p . 
Explicitly, for t — 1 we have 

e L = e Lo + e {1 - ai)La Ve aiLo da + e £-°i) L °V e^-^V e a2Lo da 

+ •••+ / e {1 - ai)Lo V . . .Ve iad - 2 - ad - l)Lo Ve ad - lL da 
J^d-i 

= e L ° + e {1 ~ ai)Lo Ve aiLo da + e {1 - ai)Lo Ve {ai - a2)Lo Ve a2L °da 

J Si J S2 

+ • • • + / e (1 "' Tl)io V . . . \Z e ( CT d-2-^-i) L o Ve ad - lLo da + ... 
J^a-i 

+ / / e (i-«n)£oy_ _ >e (a d _ 2 - CTd _ 1 )L ^ e (<Td-i-^)io^di do . d(Tn 
= e L ° + e (1 " ai)Lo 1/e ,TlLo da+ / e (1 " ol)£o ye (oi-0,,)£o ye oaio do' 

_|_..._)_ / e ( 1 -°"i) L oy e (o-i-cr 2 )Lo e (<T (i _ 1 -CT d )Loy e CT ( iL^ cr 

where each integral is well defined as a Riemann integral by the Lemma 
12321 □ 



3. Local dilations and perturbative expansions 

In this section, we tackle the task of deriving an algorithmically com- 
putable approximation to e tL . We exploit the perturbative expansion 
(I2.15P with Lq the operator obtained by freezing the highest-order co- 
efficents of L at a given, but arbitrary, point z G Mr , and dropping the 
lower-order terms (see (I3.11al) below). Then, we approximate L — Lq 
by an appropriate Taylor expansion, so that each of the terms in (12.151) 
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except the last one can be explicitly computed using commutator for- 
mulas, as discussed in Section 0] Recall that the sets of second order 
differential operators L 7 cL were introduced in Definition 12. 11 

First, using a suitable rescaling in space and time, we replace the 
problem of determining an asymptotic expansion of the kernel 

G^(x,y) :=e tL (x,y) 

of e tL by the problem of determining an asymptotic expansion of the 
kernel Q^"' z (x, y) = e 1 "" (x, y) of e LS ' z for a suitable family of operators 
L s,z parameterized by s = y/t, and by the point z G M. N . The point 
z is fixed throughout this section, but it will be allowed to vary later 
on as a function of x and y satisfying some conditions, for example 
z = (x + y)/2. For some results, we will set z — x. The family L s,z 
has limit precisely L as s — > 0. Since we will let z vary later, we shall 
sometimes write L = Lq. 

For any s > 0, we consider the action on functions of dilating x by s 
about z and t by s 2 about 0. If / : R N -> E, u : [0, oo) x 
then set 



2 c s ' z (x) 



(3.1) r z (x):=f(z + s(x-z)), 

(3.2) u s ' z {t, x) := u(s 2 t, z + s{x- z)), 

and, 

N N 

(3.3) L s ' z := "tffrWj + s b i Z ( x ) d i + s 

i,j=l i=l 

We immediately see that 

(3.4) L s ' z u s ' z = s 2 {Lu) s ' z , {d t - L s ' z )u s ' z = s 2 [{d t - L)u] s ' z 

In particular, we have the following simple lemma, which we record for 
further reference. 

Lemma 3.1. If u solves (12.101) . then u s,z solves 
(3.5) 



d t u s ' z - L s ' z u s ' z = in (0, oo) x R N 

u s,z = fs,z e c™(R N ) on {0} x R N . 



3.1. Dilations and Green's functions. We want to study the Initial 
Value Problem (13. 5p and the Green's function of its associated solution 
operator e tLS ' z . We can reduce to study the special case z = 0. 
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(3.6) 



The definition of the Green's function and Lemma 13.11 then gives 
u s >°(t,x) = [ gr°(x,y)r*(y)dy= [ Q?* \x,y)f(sy)dy 

JR N JR N 

= W gr\xj)f(y)dy. 

JRN S 



On the other hand, 

u s >°(t, x) = u(s 2 t, sx) = / G^ t (sx, y)f(y)dy, 

Jr n 

which implies 

g"*fa I) = s N gt t (sx, y) * gr°(x, y) = s N g^ t (sx, sy). 
In other words 

If we now translate to z 7^ and choose s = Vi, we obtain the 
desired correspondence between g^ and Q^ a ' , which we also record for 
further reference. 

Lemma 3.2. Assume L G L and let z be a fixed, but arbitrary, point 
in K JV . Then, for any s > 0, 

Q}{x, y) = s~ N gr(z + s-\x-z),z + s~\y - z)) 

= t~~g[ (z + t~z(x - z),z + t~*(y - z)), ifs = t~*. 

3.2. Perturbative expansion of e L "' z . Since Lemma [3.21 gives us an 
immediate procedure for obtaining the Green function g^(x, y) oidt — L 
from the Green's function g^ s ' z (x,y) of d t — L s,z , we now concentrate 
on obtaining a perturbative expansion for the latter. 

Recall that L z = L°' z = \im s ^ L s ' z . Let us write V(' z := IS* - L z . 
Then, V*' z takes the role of V in the perturbative expansion (I2.15P for 
the operator e L3 '~ , that is: 

e LS,z = e L °+ [ e ToL W°> z e TlL *dT 
+ f e T ° L °^'V lL °VfV 2i °dr +... 

+ f e ToL W°> z e TlL » . . .e Td -^ L W°' z e TdL *dT 

+ f e T ° L »V*' z e TlL ° . . . e TdL »V°' z e Td ^ LS ' z dT. 



(3.7) 
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In a sense to be made precise below, we have V{ = O(s). Consequently, 
if we let the iteration level d — > oo in (I2.15p . we obtain a formal power 
series in s. We will rigorously show in Section [5] using the exponen- 
tially weighted Sobolev spaces W£ ,p that (12.151) indeed gives rises to an 
asymptotically convergent series in s as s — ► and will derive global 
error bounds in W S)P and W* ,p for the partial sums. 

Let n G Z + be a fixed integer and consider the Taylor expansion of 
the operator L s,z up to order n in s around s = 0, 



(3-8) L s,z = £ s m L z m + K'+i 

m=0 

were is the remainder term in the expansion. Let 

T/ S > z _ B n+1 js,z 
v n+l — b ^n+l- 

The operators L z m , 1 < m < n, are given by 

1 / r/ m 

(3.9) := -r -— V 

K ' ml \ds m 

and are independent of s, while 

1 / d n+1 

(3 - 10) LT+1 := (^+1)! V^ 1 ' 

for some < a < 1, and hence it still depends on s. 

Remark 3.3. From the form of L s,z in equation ( 13. 3ft it follows that 
the operator L^, m < n, (respectively -^*'^ 1 ) has coefficients that are 
polynomials in x — z of degree at most m (respectively of degree n + 1). 
The coefficients of the polynomials themselves are bounded functions 
of z. More precisely, the coefficients of the second order derivative 
terms are of degree at most m in x — z, while the coefficients of the 
first order derivatives term are of degree at most m — 1 in x — z, and the 
coefficients of the zero order derivative term is of degree at most m — 2 
in x — z. The coefficients of these polynomials in x — z are bounded 
functions of z, together will all their derivatives, a fact that will be 
exploited later. 

The first few terms of the Taylor expansions are explicitly: 

N 

(3.11a) L z =^2 Oij(z)didj, 



TV N 



(3.11b) L1=J2 tt x ~ 2 ) ' Va ^)) didj + bi ^ z ) di 



i,j=l i=l 
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- 1 



2 

(3.11c) 

V ; AT 



+ J2((x-z)-Vb i (z))d i + c(z) 



1=1 

Since Lq has coefficients that are constant in x, from formula (12.141) 
we obtain 

(x-y) t (A )- 1 (x-y) 

(3.12) e tL 5 = - e * 

y/(Airt) N det A° 

where A := 

Furthermore V{ := L s '^ — Lq can be written as 



(3.13) V{:=^s m L 



mjz | „n+l r s,z 

S -^n+l- 

m=0 



This Taylor polynomial expansion can then be substituted into (13. 7p . 
yielding another polynomial in s. To describe each term of this poly- 
nomial and to formulate the main results in this section, we need to 
introduce some notation. Let N := {1, 2, 3, . . .} denote the set of nat- 
ural numbers (always assumed to be > 0). 

Definition 3.4. For any integers 1 < k < d + 1 and £, we shall 
denote by 21^ the set of multi-indexes a = (a±, ■ ■ ■ , a^) G N k , such 

that \a\ := Yl a j = ^- Furthermore, we denote 21^ := IJfe=i^M- ^ or 
symmetry, it will be convenient to set 21^ = {0} if £ < k, including 
when £ < 0. 

We note that, since aij > 1, the set 21^ is empty if £ < k. The 
meaning of I is that of the corresponding power of s and the meaning 
of k is that of the iteration level in the Dyson series (13. 1\ . 

Proposition 3.5. The set 21* contains 2 l ~ x elements. 

Proof. For any given, 1 < k < £, the number of elements in the set 
2lfc,£ is the number of sequences {ai,a 2 , ...ctk} of size k which add up 
to £ and is, therefore, given by ( k Z\)- Consequently, the number of 
elements in 21* is given by Yfk=\ = Efe=o (V) = D 

We are now in the position to describe the expansion 13.71 more ex- 
plicitly. We recall that d is the iteration level of the approximation and 
n is the order of the Taylor expansion. In the following definition, by 
abuse of notation, it will be convenient to write L z n+1 instead of L^+i, 
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that is, we shall omit s from the notation. We also recall that 21^ = 0, 
if i < k. This condition will be understood. 

Definition 3.6. For each multi-index a = (ati, . . . , G %lk,e, we let 

(3.14) A a , z := [ e ToL *°L z ai e^L z a2 ---L z ah e TkL *°dT, 
if 1 < k < d, and 

(3.15) A a , z : = [ e^ L o L ^ L oL z a2 ...L z ad+i e^ z dr. 

if k = d+ 1. Then, we set 

(3.16) Ai := £ A.,,, 

with the convention that A° z = e L ° . 

We observe that a uniquely determines k and £, so that our notation 
is justified. Let a = (aj) G %lk,e- We remark that if k = n + 1 or some 
atj = n + 1 (in which case stands in fact for L*'* x ), then A ayZ and 
depend on s, so we shall sometimes denote these terms by A ajSjZ 
and A^ z . 

Also, in what follows, when no confusion can arise, we will drop the 
explicit dependence on z. However, in Section \5[ z will be allowed to 
vary and we will reinstate the full notation. We also observe that each 
A z or A s z is well defined as a Riemann integral by Lemma 12.121 and by 
the following lemmas. Let us recall that (x) w = (1 + \x — w] 2 ) 1 ^ 2 . 

Lemma 3.7. The family 

{ (■''} .'!■}■ se (0,1], zeR N , j = 0,...,n + l} 

defines a bounded subset o/L. 

Proof. This is an immediate consequence of Remark 13.31 if j < n and of 
directly estimating the remainder in the Taylor series for j = n + 1. □ 

In the following Lemma, we shall use an arbitrary center for our 
weight. 

Lemma 3.8. For each given e > 0, the family 

^ e -e(z- W ) e -e(x) wL z. sG (Q 7 l] 7 zeR N , j = 0, . . . , 71 + 1} 

is a bounded subset o/L. 
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Proof. Let us assume first that w = z. We need to prove that the 
family 

{e-^L*; s G (0, 1], z G R N , j = 0, . . . , n + 1} 

is bounded in L. Indeed, this follows from Lemma 13.71 and the simple 
observation that (x)^e _e ^^ z < C, with C independent of z and j. 

To obtain the statement of the theorem, we then apply the trian- 
gle inequality to the vectors (0,x),(l,z),(l,w) G M 1+7V to conclude 

that (x — z) — (x — w) < \z — w\ < (z — w). This shows that 
ee{{x - z) - {x - w) ^ {z - w)) < L Rence the family 

^{{ x - z )-i x - w )~( z ~ w )) e ~ e ( x )z l z = e - e ( z ~ w ) e -t( x )w £ 2 j ? 

s G (0, 1], z G ~R N , j — 0, . . . , n + 1, is bounded in L, as claimed. □ 

Lemma \2 . 1 21 together with Lemma [3781 then give the following result. 

Corollary 3.9. We have A a>z G B{W°> P , W r a * e ), for any a G 2t M , 
z G M^, r, s G R, 1 < p < oo, and e > 0. Moreover, we have that 

/or a constant C q>r>p ^ e that does not depend on z. In particular, each 
A aiZ is an operator with smooth kernel A ajZ (x,y). 

Therefore, we can write 

(3.17) A a J{x)= [ A a , z (x,y)f(y)dy. 

The point of the above definition and results is to rewrite the per- 
turbative expansion (partial Dyson series) in the form 

Lemma 3.10. Denote M = (d + 1) (n + 1) . We have 

M min{£,d+l} M 

e L- = e L S + £ y: E = 

e=i k=i «e2i M t=o 

We now assume that n < d and write the perturbative expansion of 
the above Lemma as follows: 

n M 



(3.18) e LS ' z = e L o+Y,s £ Ai+ £ s £ A\ 

1=1 l=n+l 

n n 

= e L ° + J2 *'K + ," +1 E- = £ s l K + s n+1 E%, 

£=1 t=0 

where Ej 2 represents the error in the approximation and depends on 
s, whereas the terms A z , 1 < I < n do not depend on s or d, since we 
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have assumed that n < d. Since Ej* is independent of d for d > n, we 
shall eventually restrict to d = n. 

4. Commutator calculations 

The purpose of this section is to give an explicitly computable rep- 
resentation of the perturbative expansion (13.181) as 

n 

i=i 

where (x, z, d) is a differential operator with smooth coefficients that 
depend polynomialy on x — z and s, and are bounded with all deriva- 
tives in z. Both the order of the operator as well as the degree of the 
polynomial coefficients depend on the order of the Taylor expansion n, 
which also equals the iteration level d. We give an explicit characteri- 
zation of P„ and an iterative procedure to calculate it in Theorem 14.71 
The main idea is to show that each A QjZ in (I3.14p can be written as 
an explicitly computable differential operator V a acting on the distri- 
bution kernel of e L ° , and thus using (13.161) show that the perturbative 
expansion (13.181) can be rewritten in this form as well. Throughout 
this section, z is kept fixed, though arbitrary, and d will always mean 
differentiation with respect to x. 

Definition 4.1 (Spaces of Different iatial Operators). For any nonneg- 
ative integers a, b we denote by T>(a, b) the vector space of all differen- 
tiations of polynomial degree at most a and order at most b. We extend 
this definition to negative indices by defining T>(a, b) = {0} if either a 
or b is negative. By polynomial degree of A we mean the highest power 
of the polynomials appearing as coefficients in A. 

We remark that T>(0,b) consists of differential operators with con- 
stant coefficients. 

Definition 4.2 (Adjoint Representation). For any two differentiations 
A\ G T>{a\,b\) and A 2 G T>(a2,b 2 ) we define ad^^) by 

(4.1) &d Al (A 2 ) := [A U A 2 ] = A 1 A 2 -A 2 A 1 , 

as usual, and for any integer j > 1 we define ad^ (A 2 ) recursively by 

(4.2) ad^(A 2 ) ^ad^adi; 1 ^)) 

Proposition 4.3. Suppose A\ G T>{a\,bi) and A 2 G T>(a 2 ,b 2 ). Then 
for any integer k > 1, ad^^) G V(k(ai — 1) + a 2 , k(b\ — 1) + b 2 ). 
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Proof. We first notice that 

(4.3) ad Al (A 2 ) E V(ai - 1 + a 2 A - l + b 2 ). 
Next, from ( 14.21) we have 

(4.4) ad Al (A 2 ) = ad^^adyi^adA^adA^. • • )))), k - times, 

so that an application of (14. 3 p k times yields the result. □ 

Lemma 4.4. Let m,k be fixed integers > 1. Let L E P(0,2) and 
L m E T>(m, 2). Then, &d Lo (L m ) E T>{m — k, k + 2). In particular, 

(4.5) ad£ o (L m ) = 0, ifk>m. 

Proof. Applying Lemma [4.41 we see that ad^g(L^) E T>{m — k, k + 2). 
If k > m, then by definition Dim — k,k + 2) = {0} and we obtain 

(USD. □ 

Lemma 4.5. Let L$ E ^(0, 2) PI L 7 , and /et L m E V(m, 2). Tnen for 
any 9 > 0, 

L m = P m (Lo, L m ; 9, x, d)e L ° , 
where P m (9) = P m (L , L m \ 9, x, d) E V(m, m + 2) is azven fry 

m # fc 9 2 
Pm{9) := ^ — ad£ (L m ) = L m + 0[L O , L m ] + —[L , [L , L m \] + ■■■ . 

k=0 

Proof. Recall the Baker-Campbell-Hausdorff formula (see for instance 

(4.6) $(t) := e tA B - ( t k &d k A {B)/k\ J e tA = 0. 

\fc=0 / 

In general, this formula is a formal infinite series, and the equality 
$(£) = must be justified. 

Setting A = L , B = L m , we have that &d A +1 (B) = 0, by Lemma 
14.41 so the sum becomes finite, and the function $(t) is well defined as a 
bounded operator W™' p — > LP . Since $(0) = 0, to prove that $(t) = 
for all t, it is enough to show that d t §(t)f = for all / E W™ ,p . Indeed, 
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we have 
d t $(t)f = e tA ABf 

1 oo 

Y,kt k - l &d k A {B)/k\ ) e tA f- ( £f*adS(fl)/fc! ) Ae tA f 



,fe=0 / \fc=0 

oo \ / oo 



e tA 



ABf - J2 tkad A + \B)/k\ e tA f - J]t fe ad*(iL4)/fc! etA f 

\k=0 J \k=0 J 

= e tA ABf - ( t k ad k A (AB)/k\ ] e tA f 

= Ae tA Bf- A (£Yac£(£)/ifc! J e tA f = A$(t)f. 

So the continuous function u{t) := ${t)f G BP satisfies the equation 
dtu(t) — Au(t) = with initial condition w(0) = 0. By the uniqueness 
of the solutions of this equation in L p , we obtain that u(t) = 0, which 
is the desired Baker-Campbell-Hausdorff formula. 

The indicated properties of P m (9) = P m (L , L m ; 9, x, d) are obtained 
directly from Lemma POl as follows. We have ad^(-B) G T>(m — k, k + 2) 
and hence 

m nk m 

P m (0) :=J2y ad A( B ) e ^V{m- k,k + 2) C V(m,m + 2). 

k=0 ' k=0 

This completes the proof. □ 
Lemma 4.6. For a given multi-index a G 2l&^ with k < d = n, let 

V a {x,z,d):= / Y[P ai {L z ,L z a .;l-ai,x,d)da, 

where P ai (L,Q, L z a ,\ 1 — <7j, x, d) is defined in Lemma \4-5\ Then 

K,z=V a {x, z,d)e L ° 

where the product is the composition of operators and V a is a differen- 
tial operator of order 2k +£ and polynomial degree < £ = \a\ = Yli=t a *- 
More precisely, we can write 

(4.7) V a (x,z,d)=J2 E ^{z){x-zf&l 

\P\<t |7|<^+2fc 

with ag 7 G C%°(R N ) and (3 and 7 multi-indices. 

28 



Proof. The proof is a calculation based on the repeated application of 
Lemma 14.51 on A*' . We fix a £ 21^, and for simplicity we continue 
to denote P m {9) = P m (LQ, U m \ 9, x, d), when no confusion can arise. 
Then, 



A a>z = / e (1 " ffl)i »L Q1 e (ffl - ff2)i «L a2 e (ff2 " ff3)i » • • • L ak e akL «da 

= [ P Ql (l-a 1 )e (1 - <T2)L oL a2 e (CT2 - <T3)i -..L afc e CTfcL orfa 

= [ P ai (l-a 1 )P a2 (l-a 2 )e^ L °---L ak e^ L °da 
J^k 

= [ P ai (l-a 1 )P a2 (l-a 2 )---P ak (l-a k )e L °da 

k / k \ 

= [ flP ai (l-Vi)e L *°do- = I f flP^l-aJdaje^. 

The proof is complete. □ 
Finally, for I < n we set 

V\x,z,d):= £> a (x,z,0)=£; E / H p atO- 

aeVli k=l a€2t M ^ Efe i=l 

so that 



A * = E E A ^ = E E ^(^^ 

fc=l aG2t fc ^ fc=l aG2lfc,f 



EE / n^^.^il-^x.a)^ =P £ (x,^,9)e L o. 

fe=l a62l fe ,£ ^ Sfe i=l 
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A similar, but more complicated, representation holds also for A c 
and for multi-indices a e 2l n+ i^. Indeed, 

A, 



'■a.z 



e ( 1 - CT1 ) L §L ai e (<T1 - <T2)L §L a2 • ■ ■ e {an - an+l)L °L an+1 e an+lLa ' z da 



n + l 



= j P ai (l - a 1 )e {1 - a2)L ^>L a2 ■ . . e (^-n+i)Lz . . . Lan+ie °n+iL*>* da 

= / P ai (l ~ <ri) ■ ■ -P ak (l - a n+l )e^^e^ LS ' z da. 

We are now in the position to state the main result of this section. 
Below, we set V° = 1. Let us recall the error term 

(n+l) 2 (n+l) 2 n+l 

(4.8) E£- := £ s^Al = E E E 

£=n+l £=n+l k=l a£f& k: t 

introduced in Equation (13.1 8ft . (There, we introduced E^n, but such 
error term is independent of d, as long d > n, hence we can always 
assume that d = n.) 

Theorem 4.7. The perturbative expansion ( 13.181) of e LS ' z can be writ- 
ten in the form 



■ J n,ni 



e L^ =e Ll + Y^ s V(i, z , d)e L ° + s n+1 E* 

where the differential operators V e are explicitly given by Lemmas \4-5\ 
and \4-6] 

Proof. Starting with (13 . 1 8j) . we have 



n £ 

L " + EE E * e K*+s n+1 K 

1=1 k=l aGSlfc.f 



= + Sn+1 K% = + E Z > ^ + S " +lE nfn- 

e=o i=i 
This completes the proof. □ 

Recall that e L o(x,y) is explicit given in equation (13 . 1 2j) . since z is 
arbitrary, but fixed, and it agrees with the function G°(z;x,y) defined 
by equation (11.61) in the Introduction. 
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Corollary 4.8. If \a\ — £ < n, then the kernel of each operator A. ajZ 
appearing in the perturbative expansion (13.1 8p is explicitly given by: 

^ i (z,x,y)G(z;x,y), 
where the function ^ are of the form 

tf(z, x,y) = J2 «cp{z)(? ~ z) a (x - yf, 
with \a\ < I, (3 < U, a a3 e C?{R N ). 

Proof. We observe that e tL ° is a convolution operator, since z is fixed, 
therefore 

(V a {x, z, d) e L °) (x, y) = V a (x, z, d) (e L °{x, y)) . 

Then, the result follows from formula (11.61) for e L °(x,y), formula (14.71) 
for V a (x, z,d), and the fact that V is a sum of such operators as a 
varies over 21^. (See also Lemma I5TT1 in the next section, Section[5l) □ 

The method introduced in this and the previous sections to approxi- 
mate the heat kernel of L will be called the Dyson- Taylor commutator 
method. Its description is now complete. 

5. Error estimates 

In this final section, we prove all the bounds necessary to justify the 
error estimate in the asymptotic expansion of Theorem I l.ll Throughout 
this section, n will denote the order in the Taylor expansion of the 
coefficients of L, which may differ from the approximation order as 
defined in equation (13.181) . Such approximation order will be denoted 
by fi, as in the statement of Theorem \l.l\ Recall that the definition of 
the operators A a ^ z depends, in principle, on n. However, if a G 21^ and 
n is large (n > k, n> atj) the operator A a;2 no longer depends on n (in 
which case it does not depend on s either). This observation, together 
with the fact that n is fixed, justifies omitting n from the notation for 
A a)Z . Moreover, the error terms K S J Z U are independent of d, as long as 
d > v, which will always be the case, so we shall write E*'^, = E^. 
Below, we will use such error terms for v — \i and v — n, with n > \i 
appropriately chosen. 

We start from Lemma 13.101 All the terms appearing in that lemma 
are operators with smooth distribution kernels by Corollary 13.91 We 
recall that we denote by T(x, y) the distribution kernel of an operator T 
with smooth kernel (so T(x, y) is a smooth function such that Tf(x) = 
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f RN T(x,y)f(y)dy). In terms of kernels, the formula of Theorem 14.71 
takes the form 

V 

(5.1) e LS '\x,y) = e L o(x,y) + J2s i K(x,y) + s" +1 K : t(x,y) 

V 

= J2^K(x,y) + s" +1 K : t(x,y), 

where again v = fi or v = n. 

We recall that Lq is obtained from L by freezing the coefficients of 
the highest order derivatives of L at z and by discarding the lower order 
terms. 

We now substitute x = z + s~ x {x — z), y = z + s~ x {y — z), and 
z = z(x, y) in the Equation (15. 1 p above, for some function z(x, y) to be 
specified later. Lemma [3.21 and Equation (15. 1 1) then give 

(5.2) e s2L (x, y) = 8 - N e L '"(z + s~\x -z),z + s" 1 ^ - z)) 

= s£ K( z + S ~H X ~z),z + s-\y - z)) 

£=0 

+ s v+1 W£{z + s-\x -z),z + s-\y - z)), 

which is valid for any v < n, in particular for v = fi and for v = n. 

Using the definition of the approximate Green function (x, y), 
for t = s 2 , in Equation (I1.13p . we then obtain 

(5.3) e s2L {x ) y) = g [ ^\x ) y) + s^ l W^{z + s-\x-z),z + s-\y~z)). 

The error term in the approximation defined by Equation (11.131) is 
consequently given by 

(5.4) e s2L (x,y)-g [ »> z] (x,y) = s» +1 Et : l( Z + s~ 1 (x-z),z + s- 1 ( y - Z )), 

where Wjl is as in Equation (I3.18P with s = and z = z(x,y). 
We next introduce the dilated error operator 

(5.5) S [ »> z] f(x) = [ EM(z + s'^x -z),z+ s- l {y - z))f{y)dy. 

and define the approximation kernel to be the operator with kernel 
Q^' z] (x,y), so that 

(5.6) e tL _gW,z] =t (u+l)/2 s W,z]^ 

We will use the above formula only for v — // < n, where n will be 
taken large enough. 

32 



Indeed, if /i < n, then the error term can be written as, 



(5.7) E»(z + s-\x -z),z+ a-\y - z)) = 

(n+l) 2 max{£,n+l} 

E E E « < -'*- 1 Aa lJ .(2 + a- 1 (a;-z) > « + «- 1 (y-z)). 
e=/i+i k=i ct&% kt i 

(See Equation (j 1.14ft . for instance.) We will estimate Ef£* = E*'* by 
writing 

n I 

( 5 - 8 ) = E s ' _M_1 E E A «>* + sn+1_ " E n,;- 

The point of this formula is that the error term E*'^ is independent of 
n, as long as fi < n. However, splitting the error as done above will 
allow a better control on the error estimate of Theorem 11.11 In fact, 
we will show that each A e z in the first sum, which does not depend on 
s, is a pseudodifferential operator, and its contribution to the overall 
error after the parabolic rescaling will be obtained in terms of a refined 
analysis on its symbol. This analysis, in turn, leads to some refined 
estimates uniformly in s on the norm of the operator between weighted 
Sobolev Spaces. On the other hand, we will obtain only rough estimates 
on the remander term E n n , which will nevertheless be enough, due to 
the additional factor The main issue in treating the remainder 

is that some of its terms A Qi2 implicitly depend on s, a fact which 
makes it difficult to show the remainder is also a pseudodifferential 
operator, at least in the usual Hormader class. It may be possible to 
show that E n n is indeed a pseudodifferential operator employing more 
exotic symbol classes or amplitudes, but we do not need to pursue this 
point here, since we are able to prove the sharp estimates of Theorem 
11.11 in any case. 

We now proceed along these lines. The dilated error operator in- 
troduced in equation (I5.5p can be rewritten in terms of approximation 
operators 

(5.9) C s , a f(x) = s~ N [ A a , z (z + s - 1 (x-z) J z + 8 -\y - z))f(y)dy, 
as 

(n+l) 2 max{«,n+l} 

(s.io) s [ £ z] = E E E si+N '^s,a. 

t=H+l k=l a£2l fe: f 
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We therefore obtain 

n I 

(5.11) EY A = s ^~ l H H £ s , a + s n+1 ^Sl n ' z] . 

£=/j,+1 k=/i+l ae%lk,e 

To evaluate \\S /|| in a desired norm, it will then be enough to 
evaluate each operator norm ||£ S)Q! || (between suitable Sobolev spaces). 
As explained above, we shall derive a rough estimate for the terms with 
a G %l n +i,e or some «j = n + 1 (which corresponds to A a;Z depending 
on s). When A aj2 is independent of s (that is for a G 21^, k < n, 
ati < n), we shall derive some more precise estimates. We begin with 
these refined, more precise estimates. 

5.1. Precise estimates. Recall that we denote by A a)Z (x,y) the dis- 
tribution kernel of the operator A a>z since it is a smooth function. Thus, 
for a G 21^^, k < n, cti < n, A a ^ z (x,y) does not depend on s. Let us 
fix a function z(x, y), which will be specified later, and let C s>a be the 
operator with distribution kernel 

s~ N A a , z (z + s~ 1 (x- z),z + s-\y - z)), 

introduced above in Equation (15. 9p . where z = z(x,y). 

We will show below that in this range of a for a suitable choice of 
the function z, the operator £ S Q , is a pseudodifferential operator whose 
symbol is well behaved. We shall then use symbol calculus to derive 
the desired error estimates. We refer to [5T] for all relevant properties 
of pseudodifferential operators. Below, we follow the usual convention 
and set D — -d, (i — where if not specified otherwise d = d x . 

We shall need the standard seminorms p m ,a,i3 given by 

(5.12) ?w(a)= sup |(O l/3| - m <mMx,OI- 

Then the Hormander class S™ := ^q^R^ x W n ), m > — oo, is by 

definition the set of functions a : M? N — > C satisfying p m ,a,f3{a) < oo. 
The space S~°° = S~°°(j8L N x M N ) is defined by the same seminorms, 
but with m G Z arbitrary. 
We also denote by 

(5.13) Fu{x) = u{£) := [ e-^- x u(x)dx 

the usual Fourier transform of u. For any symbol a in the Hormander 
class S% := S^ (R N x R N ), we denote by a(x,D) the operator 

(5.14) a{x, D)u{x) = (2ir)~ N I e ix< a(x,£)u{£)d£, 

Jr n 
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defined for u in the Schwartz space S(R N ). We will denote by JF 2 the 
Fourier transform in the second variable of a function of two variables. 
For a E S-°° := S^(R N x R N ), the operator a(x,D) is smoothing 
with distribution kernel 

a(x,D)(x,y) = (2ir)- N f e^-^a{x, = {^a){x, x - y). 

Jr n 

Let if be a smooth function on R N x ~R N . If the integral opera- 
tor defined by K, which is smoothing, is in fact a pseudodifferen- 
tial operator a(x, D), then we can recover a from K by the formula 
(J r 2 ; 1 a)(x,y) = K(x,x- y), so 

(5.15) a(x, £)= e~ ii y K{x,x -y)dy. 

Recall next the function G(z; x) = (4vr)- JV / 2 det^^))" 1 /^-^^" 1 */ 4 
introduced in Equation (jl.6p : Then the distribution kernel of e L ° is 
given by 

(5.16) e^(x,y) = G(z;x-y), 
and we have the following result. 

Lemma 5.1. Let z eR N be a parameter and let us consider the op- 
erator T = (x — z)@d2e L ° , where (3 and 7 are multi-indices. Then the 
distribution kernel ofT is given by 

T{x,y) = {x- Z y{&lG){z-x-y). 

Proof. The Lemma follows from a direct computation. □ 

We will also need the following standard result. 

Lemma 5.2. (i) The Fourier transform in the second variable estab- 
lishes an isomorphism T 2 ■ S' 00 := S^°°(R N x R N ) -» S^ 00 . 
(ii) Multiplication defines a continuous map S^ ^ x S~°° — > S" 00 . 

(Hi) If {a s } se (o,i] is uniformly bounded in S~°° andb s (x,^) = a s (x,sl;), 
then the family {s k b s } s€ ( 1] is uniformly bounded in S±q, k>0. 

Proof. This follows from a straightforward calculation. □ 

For our main result, we require some assumptions on the dilation 
center z. 

Definition 5.3. A function z : R 2N —>■ R N will be called admissible if 

(i) z(x, x) = x, for all x G R N . 
(ii) All derivatives of z are bounded. 
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A typical example is z(x, y) = Xx + (1 — X)y, for some fixed parameter 
A. A simple application of the mean value theorem gives that (z — x) < 
C(y — x) for some constant C > 0. 

We are now ready to state and prove the main result of this subsec- 
tion. 

Theorem 5.4. Let a G %lk,e> k < n, a < n. Assume that z : 
M? N — > M. N is admissible. Then there exists a uniformly bounded family 
{a s } sG ( ,i] in S~°° such that, ifb s (x,£,) := a s (x, then 

£ s ,a = b s (x,D). 

Proof. By Lemma 14.61 we have that A a>z is a finite sum of terms of 
the form (p(z)(x — z) l3 d2e L o with <p G C£°. Let then k z (x,y) be the 
distribution kernel of a(z)(x — z) /3 d2e L ° and let 

K s (x, y) := s~ N k z (z + - z),z + - z)), z = z(x, y). 

By abuse of notation, we shall denote also by K s the integral operator 
defined by K s . It is enough then to prove our theorem for K s . Namely, 
it is enough to show that there exists a uniformly bounded family 
{a s } s£ (o,i] in S~°° such that 

K s = a s (x, sD). 

By lemma 15.11 we have that the distribution kernel of dJe L ° is of 
the form if)(z,x — y) and belong to S~°° as a function of x — y for z 
fixed. (This is consistend with the fact that for each fixed z, d2e L ° is 
a convolution operator.) More precisely ip(z,x) is ' A ^'^. 
This observation implies 

K s (x, y) = <p(z(x, y))s~ m ~ N (x - z(x, y)filj{z{x, y), - y)) =: 
ip(z)s~^~ N (x — z) l3 ip(z, — y)), z = z(x, y). 

We then let 

b s (x,£)= / e~ iy< (fi(z)s~\^~ N (x- z) p ip(z, s~ 1 y)dy, z — z(x,x-y). 
Jr n 

Next, we observe that if we change variables from y to sy, we can write 
b s (x, £) = a s (x, s£), where 

a s{x,0 = / e~ ty ' ( '(j){z)s~^{x — z) f3 tlj(z,y)dy, z = z(x,x — sy). 

We need to show that a s is a bounded family in S -00 . To this end, we 
observe that, since p G C£° and the derivatives of z are all bounded, 
q as a function of y for each x. Similarly, for each j = 
1, . . . , N, s~ 1 (xj — Zj(x, x — sy)) G Sl as a function of y for fixed x, 
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and collectively they form bounded families for s G (0, 1]. Lastly, from 
what already observed above, ip(z,y) G S~°° as a function of y for each 
fixed x. Therefore, a s G S~°° uniformly in s by Lemma \5. 21 The proof 
is complete. □ 

We now obtain the desired refined mapping property estimate by 
standard results on pseudodifferential operators. Below, t = s 2 . 

Theorem 5.5. Let a G Qi k ,e, k <n, ctj < n. Assume that z : K^xl^ 
is admissible. Then for any 1 < p < oo, any r G R, 

(5.17) ^ 2 ||£s,a||iy r >p^iy r + fe 'P < Cfe,r,pj 
for a constant Ck, r ,p independent oft G (0, 1]. 

5.2. Rough estimates. We now move to study the mapping proper- 
ties of A Q Z when either a G 2l n +i,£ or some ai — n + 1. In this case, 
the operators A a;2 depend on s also, although this dependence is not 
shown in the notation. 

The mapping properties that we establish in this subsection will al- 
low us to obtain corresponding mapping properties for the error oper- 
ator E;^, which is not immediately in the form of a pseudodifferential 
operator. Consequently, we are not able to derive bounds as those in 
Theorem 15.41 above. Nevertheless, the bounds we derive are sufficient 
to establish the sharp error estimates as t — > + in weighted Sobolev 
spaces for the overall approximation, given in Theorem ll.il This result 
is achieved by choosing judiciously an n large enough. 

As before we denote W^ p = as before, where w is the center of 
the weight (x) w = (x — w) = (w — x) used to define the exponentially 
weighted Sobolev spaces (see equation (12.41) ). We shall also write L p a = 
W® ,p . The main result of this section is the following proposition. 

Proposition 5.6. Assume that z : M. 2N —> M. N is admissible. For any 
a, any 1 < p < oo, k G Z + , r > 0, and a G 1R, 

(5.18) s \\£, Sia \\ L p^^yk,p < Cfc,p ; 

for a constant Ck, P independent of s G (0,1], of a in a bounded set, 
and independent of the center of the weight that defines the weighted 
Sobolev spaces W^ ,p . 

Proof. The proof is based on explicit kernel estimates and Riesz' lemma. 
By replacing the operator L with e a ( x ~ w ) L e - a { x - w ) ^ where w is the cen- 
ter of the weight, we can assume that a = 0, as before. 

As before, A a>z (x, y) is the smooth distribution kernel of the operator 
A a>z . For any given point v G WL N , we denote by 5@ the distribution 
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defined by S^(f) = d@f(v) (we agree that S°(f) = /(f). Then 

(5.19) d^'dl" A a , z (x,y) = (5! , (df A a , z )(S$')), 

where (, ) is the usual duality pairing. Since all the coefficients (and 
their derivatives) of L are uniformly bounded, the derivative d^A a ^ z 
will satisfy the same mapping properties as A a>z . Furthermore, for 
each multi-index j3, d^8 y G H~ q (M. N ) for q > N/2 + \/3\ and has norm 
independent of y. 

In the rest of the proof, we use the weighted Sobolev spaces intro- 
duced in (12.41) . We recall that the mapping properties between these 
spaces are uniform in term of the base point. We can therefore choose 
the weight center at x in estimating ( 15.191) . We will write = Wjf'J. 
Then 5 y G H~ q for all a G E, q > N/2 + \/3\, with 

II^ji^, := \\e a< y- x> d^ y \\ H - q < c q ^ a+ ^~ x K 

Next, we pick an e > small enough. Replacing e with e/k, where k 
is such that a G ilfc i in Corollary 13.91 yields 



— a —a — c/k 

and hence 

(5.20) \did^"A a , z (x,y)\ = \ (dP8 x >d?K,»dP"8 y )\ 

<c\\d%\\ H - q /M'KAh-^h* ji^'y*- 

— a— e' —a — a — e/fc —a 

where q > N/2 + max(\/3\, \/3'\, \/3"\). 

We will employ the bounds above to estimate 
(5.21) 

£ s , a ( x 'y) = s N A a)Z (z + s 1 (x - z),z + s 1 (y-z)), z = z(x,y). 

We first use the chain rule to conclude that, if 7 is any multi-index, 
then d2£ s , a (x,y) is a sum of terms of the form 

s- j d% <9f cf ' A a , z (z + s-'ix -z),z + s-\y - z))P, 

for appropriate multi-indices (3, /?', and j3", with P a product of factors 
of the form d a z and j < Our assumptions on z imply that p is 
bounded. Using also Equation (15.20]) . we obtain for e sufficiently small, 

(5.22) \d2C Sja {x,y)\ < Cs- N ~ hl e «(-- 1 (*-*)>-«<-- 1 (w-*)) 

< Cs -N-\y\ e -a(s-Hy- X ))/2^ z = z(x,y), 
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where the last inequality follows from (x — z) < C(y — x). From this 
inequality, we obtain after the change of variables v = s~ 1 (y — x) 



[ \d2C a , a (x,y)\dy<C a8 -M, Vie 



i N , 



and 

/ \dZC., a (x 7 y)\dx<C a8 -M, Vy£R N , 

These two estimates together with Riesz Lemma give that the map 
/ — > s^d2C s ,af is bounded from LP to L p , which is enough to establish 
the result. □ 

This proposition, and the definition of E*'^ immediately imply the 
following lemma, where as usual t = s 2 . 

Lemma 5.7. Assume that z : ~R 2N x M is admissible, then for each 
r G M, q > 0, we have 

\\£ [ t n ' Z] \\ W r, P ^ W r + q , P < C T r^/\ t G (0, T\. 

Proof. Indeed, this follows from Proposition 15.61 Equation (15.101) . and 
the continuous inclusion W r,p > L p , r > 0. For r noninteger we also 
use interpolation. □ 

We note that in the above proposition we have an additional factor of 
t~ q l 2 compared with the refined estimates of Theorem 15.51 This extra 
factor will not affect the final result, however, provided the order n of 
the Taylor expansion of L is chosen sufficiently large. 

Then, the lemma leads to the following more precise estimate for the 
error operator 8^'^. 

Theorem 5.8. Assume that z : M. 2N — > M is admissible, then we have 

W^t*' ^\\w r 'P^w r + k 'P — Crt fc//2 , t G (0, T]. 
Proof. Let us chose n+1 > fi + r and t = s 2 , as usual. Then Theorems 



5751 and 15.81 applied to Equation (15. lip give 

n I 

W^t*' ^|| W-P^W+ fe <P ^ ^ ] ^ ^ ^ ^ ^ \\£<x,z\\ W r 'P— ►W r +*,j> 

+ S n+1 -»\\£l n ' Z] \\ W r, P ^ W r + k , P < Cs~ k {l + S^ 1 "^-^) < CS 



□ 



This completes the proof of Theorem 11.11 

From (11.131) . we immediately obtain the following property on the 
principal part of the asymptotic expansion. 
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Corollary 5.9. Assume that z : M. 2N — > zs admissible. For each 
1 < p < oo, r G 1, (i > 0, and any / £ W^ p let us define 

Q [ r ] f(x):= [ g [ r ] (x,y)f(y)dy, 
Jm N 

then Qf A ! -> / in W a r ' p /or t -> 0+. 
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